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Abstract 

Annihilation processes, where the reacting particles are influenced 
by some external advective field, are one of the simplest examples of 
nonlinear statistical systems. This type of processes can be observed 
in miscellaneous chemical, biological or physical systems. In low space 
dimensions usual description by means of kinetic rate equation is not 
sufficient and the effect of density fluctuations must be taken into ac- 
count. Using perturbative renormalization group we study the influ- 
ence of random velocity fleld on the kinetics of single-species annihila- 
tion reaction at and below its critical dimension dc = 2. The advecting 
velocity field is modelled by the self-similar in space Gaussian variable 
finite correlated in time (Antonov-Kraichnan model). Effect of the 
compressibility of velocity field is taken into account and the model 
is analyzed near its critical dimension by means of three-parameter 
expansion in e, A and 77. Here e is the deviation from the Kolmogorov 
scaling, A is the deviation from the (critical) space dimension 2 and r] 
is the deviation from the parabolic dispersion law. Depending on the 
value of these exponents and the value of compressiblity parameter a, 
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the studied model can exhibit various asymptotic (long-time) regimes 
corresponding to the infrared (IR) fixed points of the renormalization 
group. The possible regimes are summarized and the decay rates for 
the mean particle number are calculated in the leading order of the 
perturbation theory. 

1 Introduction 

Variety of chemical reactions occur in fluid environment and represents very 
important subject in various chemical, biological or physical systems [III21I3]- 
In these processes the reacting particles are affected not only by the diffusion 
motion but also by the external fluid flow. The usual approach to such 
reactive flows is based on some combination of reaction-transport equations. 

In this work we will concentrate on the study of the annihilation reac- 
tion A + A ^ 0, which is the paradigmatic model for the reaction-diffusion 
processes. For this type of reaction in low space dimensions the usual de- 
scription by means of kinetic rate equation is not sufficient and the effect 
of density fluctuations must be taken into the account [1]. It can be shown 
that the upper critical dimension for this process, even in the absence of 
advective flow, is two due to the density fluctuations. A renormalization 
group treatment was successfully applied for different choices of reactive flow, 
e.g. time-independent random drift [3], velocity field described by stochastic 
Navier-Stokes equation [6] or short-range correlated potential disorder [7]. 
Often large differences are observed in relation to the experiment and their 
origin can be caused by the presence of large-scale anisotropies, effect of com- 
pressibility or parity violation. There are important differences [SI E] between 
advection of scalar quantities like density on one hand and like tracer (tem- 
perature, concentration) on the other by compressible versus incompressible 
flow. It was shown that compressibility could lead to the slowing of trans- 
port process for scalar admixture and also to the enhancement of intermittent 
phenomena. These effects might be understood as a result of inhibition of 
separation between particle trajectories and therefore we expect that react- 
ing particles would spent effectively more time in the mutual vicinity than 
in the incompressible case. Hence in this case we expect the faster decay 
rate than for the incompressible case. Hence it would be desirable to sup- 
port such naive picture in more quantitative manner. In order to do that 
we will apply Antonov-Kraichnan model pUj [TT] for describing advection of 
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reactive A particles and we will present renormalization group (RG) study in 
the vicinity of its critical dimension dc = 2. In the one-loop order all relevant 
physical quantities are calculated and contrary to the incompressible case it 
is found out, that already in the one- loop approximation fluctuations of the 
velocity field affect the renormalization of the rate constant. 

After brief description of the model in Sec. [H results of RG calculations 
are presented in Sec. O In Sec. S] possible large-scale regimes are listed and 
their physical interpretation is concluded in Sec. |5l 

2 Field-theoretic model 

Field-theoretic action for the annihilation reaction process 



can be obtained in the straightforward way [12] employing Doi approach [13] . 
It can be written in the following standard form 

+ no y dx^t(x,o), (1) 

where Dq is the diffusion constant of reacting particles, and because of dimen- 
sional reasons we have rewritten product XqDq instead of the rate constant 
Kq. Last term in the action stands for the initial conditions, which are 
traditionally chosen in the form of Poisson distribution. 

In compressible velocity field, there are two types of diffusion-advection 
problems: advection of a density field and advection of a tracer field 
Here, the case of advection of a density field will be analyzed [13]. In the 
action ([1]) this corresponds to the replacement ip'^dfip — )■ ip'^ [dtip + (V.v^)], 
where v = v(x, t) is the advecting velocity field. According to [11] let us as- 
sume that V is a random Gaussian variable with zero mean and the correlator 
(in the frequency-momentum representation) 

„ r)3t.2-2A-2e-2»y 

(^o.^o,)o = ^;/(^^^^p_,), [P..(k) + «Q.,(k)], (2) 

where go is the coupling constant, the exponents e,A and r] play the role of 
small expansion parameters. They could be regarded as an analog of the 
expansion parameter e = 4 — d used in the theory of critical phenomena. 
However, in this paper e should be understood as deviation of exponent of 
the power law from that of the Kolmogorov scaling [TH], whereas A is defined 
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as the deviation from the space dimension two via relation d = 2 + 2 A, and 
the exponent t] is related to the reciprocal of the correlation time at the wave 
number k. The parameter serves for labeling of the fixed points and it can 
be interpreted as a ratio of velocity correlation time and the scalar turnover 
time [T7]. In ([2]) besides the standard (incompressible) transverse projection 
operator -Pij(k) = 5,^ — kikj/k"^ the longitudinal projector (5jj(k) = kikj/k'^ 
has been introduced. The positive (necessary for positive definiteness of the 
correlator (vv)) parameter a represents the degree of compressibility. The 
incompressible case is obtained by the setting a = 0. 

The Antonov-Kraichnan model for the advection field v contains two cases 
of special interest: 

(a) in the limit Uq oo, = Qq/uI = const we get the 'the rapid-change 
model' D^{uj, k) — )■ g'^Dok"'^"'^^''^'^^'^ ^ which is characterized by the white- 
in-time nature of the velocity correlator. 

(b) limit Mo — 0,(7q = Qq/uq = const corresponds to the case of a frozen 
velocity field Dy{u,k) gQD'^7i6{uj)k'^^~'^^, when the velocity field is 
quenched (time-independent). 

The averaging procedure with respect to the velocity field v(x) may be per- 
formed with the aid of the following action functional 

82 = -^ J dtdx J dt'dx v(t, x)D;\t -t',x- xXt', x), (3) 

where is the inverse correlator ([2]) (in the sense of the Fourier transform). 
The expectation value of any relevant physical observable may be calculated 
using the complete weight functional yV{ip'^ ,ip,^) = e^'^'^^^, where 5*1 and S2 
are the action functionals ([H) and ([3]). 

3 UV renormalization 

The inclusion of longitudinal part Q into the correlator for velocity field 
does not affect the renormalization group analysis developed for the such 
model [18]. Therefore we just mention main steps of theoretical description 
and deviations from it caused by compressibility violation. All canonical di- 
mensions of fields and parameters are listed in Tab. [H The only difference 
with the incompressible case is that now the velocity field has to be renor- 
malized [TT] . Following [21 HH HB] it is easy to prove that the model under 
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Table 1: Canonical dimensions for the parameters and fields of the model 



consideration is multiplicatively renormalizable and can be made UV finite 
by the following renormalization prescription 



with the additional constraints between them 



^g^D ~ ZuZf) — 1, ZgZ^ — Z'^, 



(4) 
(5) 



which are the consequences of the absence of renormalization of non-local 
term (|2]). Non-local character of the velocity correlator is caused by the 
nontrivial correlations in momentum and frequency scales. The total renor- 
malized action can be written as 



dt 



^ - ^WZdV^^ + ^^Z„{V.v^)] 



dtdx / dt'dx' 



+ 



'0^1 

+ no y dx ^"^(x,0). 

The perturbative calculation of the renormalization constants in dimensional 
regularization with the use minimal subtraction (MS) scheme is straightfor- 
ward [in] . We restrict ourselves to the first order in perturbation theory and 
this approximation already contains first nontrivial effect of the compress- 
ibility. It can be seen from the perturbation expansion of the one-particle 
irreducible function T^^^2 (known as interaction vertex) 



ADXZxfi' 




(6) 
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where higher order terms in couphng constants were neglected. The first 
graph on the r.h.s. of physically represents process of density fluctua- 
tions that account for annihilating of two particles. It is easy to see that 
transverse character (equivalent statement to the incompressibility condition 
for velocity field v) of propagator (vv) leads to UV divergent contribution of 
the second graph in (E]), whereas for the models with incompressible velocity 
field V it leads to UV convergent contribution. Physically this graph can be 
interpreted as an attracting process that brings together particles into a sink 
of compression, which can lead to effective increase of the reaction rate as 
will be pointed later. In the MS scheme the renormalization constants in 
one-loop calculation obtain the following form 



Zd 
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A ag 
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(7) 
(8) 

(9) 



and from relations (j5]) constants Zu and Zg can be calculated. Limiting case 
a = agrees with the results (to the one-loop precision) for incompressible 



case [18] and case g = leads to the presence of only density fluctuations [1]. 

From the relations and the beta functions f3g, and f3x are ob- 
tained via the standard definition (]g = fid^g\o (the subscript "0" refers to 
partial derivatives taken at fixed values of the bare (unrenormalized) param- 
eters) 

(3g = g[-2e-7] + 3-fD~2^.], A. = ^[-77+71)], /3a = A[2A-7a+7d], (10) 

where the anomalous dimensions '-/p are defined as [19] 

If = /i^M In Zf\o = iPgdg + P^du + Pxdx) In Zp. (11) 

It was conjectured [HI [T71 120] , that up to the two-loop calculations there is no 
direct influence of the parameter rj on the anomalous dimensions. Hence one 
can use in the actual calculations of the Feynman graphs different values of rj. 
In our calculations the simplest choice 77 = was apphed. Finally substituting 
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into the definition (1111) the anomalous dimensions are obtained 
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4 IR stable regimes 

We are interested in the IR asymptotics of small momentum p and fre- 
quencies u of the renormalized functions or, equivalently, large relative dis- 
tances and time differences in the {t, x) representation. Such a behavior 
is governed by the IR-stable fixed point g* = {gl,u* , \*), which are de- 
termined as zeroes of the (3 functions (ITUI) : f3{g*) = . It is said, that 
the fixed point g* is IR stable, if real parts of all eigenvalues of the matrix 
flij = d(3i/dgj\g=g*; i,j e {g,u, A} are strictly positive. 

The simplest way to find the average number density n{t) = (ipit)) is 
to calculate it from the stationarity condition of the functional Legendre 
transform [2T] of the generating functional obtained by replacing the un- 
renormalized action by the renormalized one in the weight functional [B]. 
This is a convenient way to avoid any summing procedures used [1] to take 
into account the higher-order terms in tiq. For a spatially homogenous solu- 
tion this leads to the rate equation with the initial condition n(0) = no for 
the average number density n(t) = {'ipit)) 



where no is the initial number density. Since the fields ip and -ip^ are not 
renormalized, the Callan-Symanzik equation for the mean particle number is 
easily obtained by the standard procedure [6l 122] 



Solving it by the means of characteristics it can be shown (details will be 
published elsewhere [23]) that the value of the decay exponent (3 defined 
through the asymptotic relation : n{t) ~ is given by the expression 



nit) 
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/3 = 1 + 
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Note that in contrast to the previous studies [HI EB] here we have to deal 
only with the three-charges' {g, u, A} theory. Passive character (no backward 
influence on the advecting fleld) of the reacting particle manifests itself also 
in the relations dx/3g = dx(3u = resulting from (|9]) and f lTOj) . These relations 
greatly simplify calculation of the eigenvalues of Uij matrix. 

Detailed analysis of flxed point structure reveals that studied system can 
exhibit one of ten possible IR regimes listed below. First let us consider 
the "rapid-change model" {u — )■ oo) . Introducing convenient variables w = 
l/u,g' = g/v?., the corresponding /3 functions can be written in the form 



where anomalous dimensions are 



w[r]-iD\, /3a = A[2A - 7a + 7d], 
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The "rapid-change model" corresponds to the flxed point with the value 
w* = 0. In this case four stable IR flxed points can be realized: 



FP 1: 



FP 2: 



FP 3: 



FP 4: 
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For the analysis of the regime m — (quenched velocity field) we introduce 
the new variable g" = g/u. Hence the corresponding (3 functions have the 
form 

= /[-2e + - 27J , /3„ = n[-r/ + 7z,], /3a = A[2A - 7^ + 7d] • 
and anomalous dimensions are given as 
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The quenched regime corresponds to the m* = and also in this case there 
are four possible IR stable fixed points: 
J I* 



FP 5: 



FP 6: 



FP 7: 



FP 8: 



g- =0, A* = 0; 
VLi = -2e, ^2 = -r], ^3 = 2A 
13 = 1; 

g"* = 0, A* = -47rA; 

ni = -2e, ^2 = -7], ^3 = -2A; 

/3 = 1 + A; 

g"* = 8ne, A* = 0; 

Q^=2e, ^2 = -^7 + e(l -a), ^3 = 2A + e(l - 2a); 

2 - (1 - 2a)6 
^ 2-(l-a)e' 

/* = 87re, A* = 27r[-2A + e(2a-l)]; 

= 2e, ^9 = -r/ + (1 - a)e, ^3 = -2A + (2a - l)e; 



2- (1 -a)e 

The nontrivial case occurs when no special choice for parameter u is con- 
sidered, i.e. let's consider u finite and non-zero. Solving equations ffTOj) for 
u ^ 0, g ^ the following values for the coordinates of the fixed point are 
obtained 



2e-r] 
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a{ri - 2e) 



(12) 



87rn*(l + n*) 1 + a' (l + a)(7/-e)' 

The two possible regimes are distinguished by the value of the coordinate A*. 
Fixed point with zero value is given by 

2-ri + a{l + e-7]) 



FP 9: A* = 0, P 



(l + a)(2-r/) 



and is stable in the region 



(l-a)6<r/<e, A + ^— — ^ > — - . 13 

2(1 + aj 1 + a 

The fixed point with non-zero value of A* is given by 

A 27rae 1 + 2a ^ 2 + 2A 

FP 10: A* = -47rA + - 27rr/-^, /3 = , 

l + a 1 + a 2 — ?7 

and is stable in the region 

2(1 + a) 1 + a 

5 Conclusions 



Fixed points 1 and 5 corresponds to the non-interacting (mean-field or 
Gaussian) theory and thus their predictions should agree with the ones of 
rate equation approach, which is indeed the case [1]. 

From the fixed points' structure some physical consequences can be de- 
duced. First we see, that compressibility has direct infiuence on the value 
of decay exponent (see FP 3, FP 7-9). For some regimes (FP 3, FP 7) 
it could lead to the enhancement (for both of them /3 > 1) of the reac- 
tion process compared to the corresponding regimes for incompressible case 
[6], [18] . As was already pointed this fact can be explained by the presence 
of compressible sinks into which particles are attracted (see also Sec. 4 in 
[2l]). However we also observe that when both density fiuctuations and com- 
pressibility are relevant (FP 4, FP 10), the density tends to suppress the 
infiuence of compressibility. 

The "real problem" corresponds to the choice e = ?7 = 4/3, which leads to 
the famous Kolmogorov "five-thirds law" [16] for the spatial velocity statis- 
tics. It is easy to see, that this regime can be realized either by the fixed 
point FP 9 or FP 10 depending on the value of parameter A, or equiva- 
lently on the space dimension d. From ( fT3l) and (fT^ we see that there is a 
"critical" value Ac = —2/(3 -|- 3a) for the parameter A, above which FP 9 
is stable, whereas below it FP 10 is stable. Because a is positive quantity, 
a parameter Ac is negative. Therefore in the vicinity of the space dimension 
two (A = 0) regime represented by the fixed point FP 9 should be realized 
with the decay exponent /3 = (1 + 3a;)/(l + a). Thus we can conclude, that 
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compressibility has a profound effect on the large-scale asymptotic behavior 
of the annihilation process and leads to the enhancement of it near its critical 
dimension dc = 2. On the other hand density fluctuations leads again to the 
suppression of compressibility (as can be seen by direct numerical inspection 
of exponent (3 for FP 10). 
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